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(57) Abstract 
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METHOD AND SYSTEM FOR MAPPING TRAFFIC CONGESTION 

FTF.I.D OF TMF. TNVENTION 
Thisinventionrelatesgenerallytoamethodandsystemformappm^ 
and in particular to a method for improving the accuracy of said mapping when a relatively 
5 small percentage of vehicles are used as traffic probes. 

ArTCOROI TNP "F TT4P. TNVENTION 
Traffic congestion is an increasingly serious problem in cities. 
One way to identify and map such congestion in real time (the first step to relieving 
it) is to identify and map the positions of vehicles that are stopped or moving slowly. Such 
10 systems are often referred to traffic control and car navigation in the field of Intelligent 

Transport Systems (ITS). 

PCT publication WO 96/14586, published 17 May 1996, the disclosure of which is 
incorporated herein by reference, describes, inter alia, a system for mapping of vehicles in 
congestion. 

15 In one embodiment described in the above publication, a central station broadcasts a 

call to thevehicleswhichrequests those vehicles which are stopped or which have an average 
velocity below a given value to broadcast a signal indicative of their position. Such signals 
are broadcast in slots, each of which represent one bit (yes or no) which relates to a position. 
Preferably, only one logical 
20 to define the related position. Such signals are then used to generate a map of those regions 
for which traffic is delayed or otherwise moving slowly. 

Preferably, an additional call is sent to the vehicles requesting transmission of 
indication signals which locate the slow moving or delayed vehicles at a higher resolution 
than that of the first call. Further calls may be made to allow for transmission of additional 
25 information on the status of the vehicles and/or to provide further characterization of the 
delays. 

Fig. 1 shows an initial map generated by such a method, wherein the area represented 
by a pixel (slot) may, for example, be of the order of 250 to 1000 meters square. 
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In a preferred embodiment of the invention described, the system then determines, 
based, inter alia, on the extent of the various contiguous areas which shows positive 
responses, a smaller area or areas for further study. Preferably, the system broadcasts a further 
query requesting those vehicles within the smaller area that have at least a given delay (which 

5 may be the same as or different from that used in the first query) to broadcast in slots, each 
representingaposition,usingafiner resolution, for example,100 to 250 meters square. Based 
on the responses to this query a second map such as that shown in Fig. 2 is generated. As can 
be seen from Fig. 2, various branches of a road network radiating from an intersection, 
designated as A-F in Fig. 2, can be identified. To improve the usefulness of the display, a 

10 background map, such as a road map, may be displayed underlying the displays of any of 

Figs. 1, 2 or 4 (described infra). 

In the event that additional information relating to the delay is desired, further queries 
can be made. For example, vehicles which are traveling toward the intersection can be 
requested to broadcast in a slot which corresponds to the slot they are in and to their velocity 
15 toward the intersection. This allows for generation of the graph shown in the lower portion 
of Fig. 3. Additional slots may be used for the acquisition of other information regarding the 
responding vehicles. Such information may also be graphed as shown in the upper portion of 
Fig. 3. 

Alternatively or additionally, a map which shows the average velocity of the vehicles 
20 toward the intersection as a function of the position can be generated. Such a map is shown 
in Fig. 4. To acquire the information needed for generating such a map, a number of queries 
maybemade, eachrequesting an to 

a given average velocity toward the intersection. The responding vehicles would broadcast 
their indication signals in slots corresponding to their position. In the map of Fig. 4, the 
25 velocity for a given pixel is determined, for example, as the average velocity of the reporting 
slots for that position. In a display of the map of Fig. 4, the velocity or delay toward the 
intersection can, for example, be displayed as a gray scale value or as a color, with for 
example red being the highest velocity or delay andbluebeing aminimum displayed velocity 
or delay. 

30 Fig. 5, is a generalized block diagram for a system useful for performing the ITS 

function described above (and which is also useful for the method of the present invention). 
Fig. 5 shows a base station or control center 9 1 having a control center transmitter 79 which 
broadcasts queries and optionally other signals to vehicles on command from a control 
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computer 80. A remote vehicle 85 (only one vehicle is shown for simplicity) receives the 
query at a vehicle receiver 84 and transmits commands to a microprocessor 86, based on the 
queries it receives from the control center. 

Microprocessor 86 also receives information regarding the status of the vehicle from 
5 one or more information generators and sensors indicated by reference numeral 88. This 
information may be sent by the sensors on a regular basis or may be sent on command from 
the microprocessor. 

Microprocessor 86 is then operative to command vehicle transmitter 90 to transmit 
indication signals (or if required, information bearing signals) in a suitable slot in accordance 
1 0 with the information received by microprocessor 86. 

The indication (or other) signals are received by a control center receiver 92 and 
processed by receiver 92 and computer 80. While the operation and construction of the 
apparatus designated by reference numerals 82, 84, 86 and 90 is straightforward and needs no 
further explanation, the operation of receiver 92 is usefully expanded upon with reference to 
15 Fig. 6. 

Generally speaking, the RF signals transmittedby the vehicle may be at any frequency 
slot. It is to be expected that there will a certain amount of frequency diversity caused by the 
imperfect accuracy and stability of the vehicle transmitters 90. The slots are wide enough to 
accommodate this diversity. 
20 Furthermore, often the system utilizes very large numbers of vehicles. If too many of 

these vehicles (in some particular situation) transmit in the same slot, then the total power 
transmitted may exceed authorized ERP or dynamic range restrictions. To overcome this 
problem longer, lower power, pulses may be used for indication signals. Furthermore, if a 
single receiver is used for receiving signals for all of the slots, intermodulation effects may 
25 cause spurious signals to appear in slots for which no actual signals have been received. 

These problems aswell as near-end to far-end transmission problems are substantially 
solved by the system shown in Fig. 6 and by certain constraints placed on the system which 
are not shown in Fig. 6. The problems and constraints but are described in the above 
referenced PCT publication, which should be consulted for a more complete exposition of the 
30 method and apparatus shown in Figs. 1-6. 

Fig. 6 shows a receiver system corresponding generally to reference number 92 and 
to a portion of computer 80 of Fig. 5. While the system of Fig. 6 is suitable for the prior art 
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US system ofthe PCT publication, it is also suitable for use with the ITS system ofthe 
present invention. 

An antenna 94 (or an array of antennas) receives signals from a plurality of vehicles 
simultaneously and passes them to a receiver and (optionally) AGC 96. Receiver and AGC 
5 96, which may be of conventional design, downconverts the received signals from RF to IF 
frequencies. The threshold levels ofthe detection process may be dependent on the AGC 
process. Alternatively, the system is operated in a closed loop mode in which the power 
radiated by the vehicles is dependent on the power received by the base station. 

The IF signal is digitized by art A/D system 98 and further down converted by a 
10 downconverter 100 to base band. It should be understood that this receiver/downconverter 
system does not demodulate the incoming signals, but only downconverts the RF so that the 
same relative frequency differences ofthe signals is present at the output of converter 100 as 
in the incoming signals, except that the absolute frequency has been reduced to a low 
frequency from the RF frequency ofthe transmitted signal. At these lower frequencies digital 
1 5 systems can be used to analyze and detect the signals. 

The low frequency band signals are fed to a series of correlation filters 102 

(correlation-typerec^ 

correlation time of the correlation filter. Preferably, the frequency bandwidths of adjacent 
receivers 102 overlap so that the entire bandwidth of each ofthe slots is covered by one set 
20 of receivers 102. The output of each of the receivers is compared to a threshold 104 to 
determine if a signal is present at the frequency ofthe respective receiver 1 02 and the outputs 
of all of threshold detectors for a given slot are OR gated (or the best signal is selected) to 
determine if any signal is present in the slot. 

In an alternative preferred embodiment ofthe embodiment disclosed, the strongest 
25 output ofthe set of correlation receivers is chosen for comparison with a threshold, with or 
without post-detection integration. 

Use of a plurality of overlapping narrow band receivers in this manner also reduces 
the extent of side lobes ofthe detection process outside the band ofthe slot. This allows for 
closer frequency spacing of the slots since interference between slots having adjacent 

30 frequencies is reduced. 

One set of receivers 102, threshold detectors 104 and an OR gate is provided for each 
slot and is referred to herein as a slot detector unit. Slot detector units for all ofthe slots feed 
adataprocessorlOSwhich, together with computer 80 processes the data as described above. 
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When large numbers of vehicles are used in the system and intermodulation becomes a 
problem (or if AGC is used, and low level signals are lost), it may be necessary to provide a 
plurality of front end portions of receiver 92 (the front end being defined as receiver 96, 
converter 98 and converter 100), where each front end receives signals from only a portion 
5 of the entire frequency band including one or many of the slots. The function of correlation 
receivers 102 may also be implemented, for example, using set of DFT's or an FFT (for CW 
signals), matched filters or other correlation receiver methods or other optimum receiver 
methods, depending on the transmitted signals. Othermethods such as energy detectors (e.g., 
radiometers) with or without tracking may also be used, however, they will give less optimal 
10 results, because of practical limitations on input band-pass filter designs. 

It should be understood that using a plurality of correlation receivers for the same slot 
may increase the false alarm probability and hence the threshold for positive detection may 
be adjusted to provide a desired low false alarm probability. 

The system may also be provided with a display 1 10 for displaying the data, and with 
15 a user interface 1 12 which is used by an operator to control both the operation of the system. 
The user interface also preferably controls the display and the memory to allow for the 
operator to review the maps previously generated or to generated new displays based on 
information previously received. 

This system works well. However, there is a need for improved accuracy of mapping 
20 and/or utilizing a relatively small percentage of participating vehicles. 

ST TMMARV HT? THF INVENTION 
The present invention provides a system and method for mapping parameters of a 
traffic congestion, e.g., a road congestion, relative to a focus. Mapping of the road 
congestion may include determination of an average length of the road congestion over a time 
25 interval, motion rate in the road congestion and arrival rate to the road congestion. These 
parameters, in turn, may be used to determine an expected delay in traveling throughout the 
road congestion as well as trends (i.e., changes with time) in the road congestion. 

The mapping is performed relative to a mapping focus, typically the front end of a road 
congestion. Themapping focus is preferably identified using the system and method described 
30 in the above-mentioned PCT Publication WO 96/14586, or it may be identified using any 
other suitable method known in the art, for example, by simple polling of predesignated target 
vehicles. Alternatively, the mapping focus may be provided from an external source, for 
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example, based on reports regarding a problematic intersection or suspect intersections that 
are to be continuously monitored. 

In an embodiment of the present invention, the mapping system constructs snapshots 
of mapping samples received from a small percentage of predesignated probes, e.g., a small 
5 percentage of vehicles equipped with an appropriate receiver and transmitter. The mapping 
samples are preferably received in response to predefined broadcast queries sent from the 
mapping system. In an embodiment of the invention, the determination ofthe average length 
of a road congestion may be based on a direct approach, obviating the need to estimate 
discrete lengths ofthe road congestion, in dynamic conditions that may include variations in 
10 the arrival rate of vehicles to the road congestion and the departure rate of vehicles from the 
congestionovertime. In preferred embodiments ofthe invention, the average motion 

ratewitmnmeroadcongestionmaybedeterminedwithoutmeneedtomcre^ 
of the radio spectrum that is used by the probe vehicles. The determination of motion rate m 
addition to the length ofthe congestion enables to estimate the expected time delay for a 
15 vemclematisabouttoentermeroadcongestionThememodofdetermimngmo 

isolatereportervehicle^^ 

invention, e.g, to increase the accuracy at which the length of the road congestion may be 
determined and/or to provide information about trends intheroad congestion, i.e., expansion 
or contraction ofthe congestion. Such isolated mapping enables the system, for example, 
20 to concatenate non-overlapping segments in the mapping samples and, thereby, to estimate 
theaveragearrivalratetome^ 

an estimated departure rate to provide trends in the average length over time, whereby a 
preferred path chosen by a vehicle may be selected based on a current time delay as well as 

on the trend in the road congestion. 
25 The concatenation of non-overlapping segments of mapping samples, in accordance 

with a preferred embodiment of the invention, may also be useful for estimating the 
percentage probe vehicles within the road congestion. Based on this estimation, in 
conjunctionwithacalculationofthe estimated arrival rate and the estimated motion rate, the 
average length may be determined more accurately. This means that the number of mapping 
30 samples canbe optimized to provide an accurate determination ofthe average length ofthe 
road congestion based on the parameters described above. Pre-stored data which may be 
generated by computer simulation of different road congestion conditions may be used in 
determining the optimum number of mapping samples for average length determinahon. 
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In accordance with the present invention, as described herein, concatenated mapping 
samples may be used to estimate the arrival rate and the percentage of probes. In traffic 
situations where two or more road congestions are correlated, several such concatenations 
fromseverdrifferentroad congestion™ 
5 example, to estimate the percentage of probes based on Maximum Likelihood estimator for 
Binomial distribution, the concatenation of more than one concatenated mapping samples 
from several correlated roads may be used by a statistical estimator. This can be used to 
improve estimates from short concatenated sample at an early stage of mapping a road 
congestion. 

10 The motion rate within the road congestion, which may be detected based on two 

mapping samples, may also ne used for determining a minimum required rate for taking 
snapshots of mapping samples according to a desired accuracy in determining the average 
congestion length. In embodiments of the present invention, the level of accuracy in 
determining average length based on motion rate may be estimated by computer simulation 
15 and provided as pres-stored data to determine an appropriate mapping sample rate. As 
mentioned above, motion rate can be detected by two mapping samples. At an initial stage 
of a mapping process, when the average arrival rate of vehicles to the road congestion and the 
probabilityofanarrivingvehiclebeingaprobecannotbe correctly calculated, prior statistical 
data may be used to initiate the estimation process. Refinement of these initial values may be 
20 performed during the sampling process by constructing concatenated segments of non- 
overlapping mapping sample segments and determining the average arrival rate as well as the 
percentage of probes, thereby enabling to determine the probability of an arriving vehicle 
being a probe. A similar approach may be used for determining the number of mapping 
samples according to the pre-stored data. 
25 The pre-stored data may be based on computer simulation to provide minimum error 

inthedeterminationofthe average length or a modified average length. The modified average 
length may take into account predetermined parameters, e.g., giving more weight to later 
mapping samples than to earUer mapping samples or any other desired criteriathat may result 
inamore accurate estimation process. As traffic condition are being mapped, statistical data 
30 is collected relating to average arrival rates and distribution of probe vehicles, whereby the 
system converges to realistic values at relatively early stages of the mapping, even before one 
would expect to have sufficient mapping samples to estimate these parameters. 

In case of traffic light control, the sampling rate may be adjusted in accordance with 
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the rate of change of the traffic lights, e.g., the timing of the green light activations. The 
timing of light changes may be provided by probe reports according to their reaction time to 
green light setting calibrated to distance from the traffic light. According to this embodiment, 
the times may be provided by a report from a probe which has been isolated for the purpose 
5 of motion rate estimation and other estimations, as described above. It should be noted that 
me average road congestion len^^^ 

rate in each mapping sample is substantially equivalent to the average arrival rate. 

When the average departure rate is not equal to the average arrival rate, the departure 
rate may be artificially adjusted to increase or decrease the length of the mapping samples, 
10 thereby to adapt the average departure rate to the average arrival rate. This may assist in 
determining the length of road congestion. Once the road congestion length is determined, 
based onthe artificial adjustment, ^adjustment stagemay be applied to compensate for the 
artificial adjustment. The compensation maybe based on a new weighted average which takes 
into account a trend in the length of road congestion. At any given time, the average length 
15 determination may be based on the latest mapping samples according to the number of 
mapping samples that will best determine the average length of the road congestion. 
Successive average length values may fed through an appropriate filter,as is know in the art, 
to remove large, random changes in value. 

The present invention is comprised in a number of improvements on the prior art 
20 system which improve the position related accuracy of the system. 

As in the prior art system described above, preferred embodiments of the present 
invention may utilize the position related data transmission system of the above referenced 
PCT publication. In addition, the present invention may utilize the general structure of the 
transmitter andreceiver as described m mat publication and in the Background of the present 
25 invention.Itshouldbenotedmat,becausemepresentmventionmayutilizeacommunication 
platform and related technology similar to those described in the above mentioned 
publications,manyfeaturesofthemethods, devices and systems described in that publication 

are also applicable to the present invention. 

Accordingtosomeaspectsofsomepreferredembodimentsofthe invention, mapping 

30 of congestion is based on identification of the starting point of traffic congestion and a 
determination of a distance of vehicles from a congestion start point or focus. The length of 
the congestion is estimated from the distance of the vehicle farthest from the congestion 
whose velocity is below a given velocity, preferably for some minimal time period. 
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Preferably, the vehicle positions are not determined for individual vehicles. Rather the 
vehicle report according to their positions, that correspond to a predetermined sub-area, if 
they are stopped or if their velocity is below some value. 

Preferably, vehicle positions over a time period are combined to form a congestion 
5 map. Preferably, the positions that are combined are determined at the same position 
resolution. Alternatively, they do not. 

According to an aspect of some preferred embodiment of the invention, the position 
of a vehicle is reported based on a distance to a known focus of a congestion. 

In a preferred embodiment of the invention, the location of a potential congestion is 
10 determined by vehicles that are stopped or moving slowly reporting their positions at a low 
resolution, for example using a rectangular grid for two dimensional mapping. Once a 
potential congestion is identified, the position of the vehicles is reported based on their 
distance from a focus of congestion. 

There is thus provided, in accordance with a preferred embodiment of the invention, 
1 5 a method of estimating the position, in an ITS system, of the length of congestion at a focus 
of a slowdown, the method comprising: 

determining the positions of one or more vehicles farthest from the focus as a function 

of time; and 

estimating the length of the congestion based on the function. 
20 Preferably, thepositionisestimat^ 

Preferably, the position is estimated as the position of a vehicle furthest from the focus 
during a given preceding time period. 

There is further provided, in accordance with a preferred embodiment of the invention 
a method of improving the reliability of an ITS system, comprising: 
25 determining the position of a plurality of vehicles; 

determining an indication of a traffic stoppage if more than one vehicle is stopped 

along a line of vehicles. 

Throughout this disclosure, where applicable, the terms and phrases listed below may 

be denned as follows: 
30 Ma pping focus : 

A position in a mapped road that defines the front end of the mapping range towards 
traffic moves usually refers to the front end of a road congestion. 
Probe : 
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A vehicle equipped with a transmitter connected to a computer both comprising an 
intelligent transmitter wherein the computer is provided with timing and positioning 
peripherals that according to a predetenruned procedure can identify congested conditions^ 
motion cycles parameters in a congested road, preferably equipped also with a receiver that 
5 enables a mapping system to control the activity of the reports preferably including levels of 
congestion to be experienced by the probe fore a report, resolution of position report, actual 
report time of a characteristic value of its position, disabling transmission of probes that are 
closer than a certain position to the mapping focus and re-enabling the transmission, and 
according to a predetermined protocol reports will preferably include, but not limited to, one 

10 or more of the following: 

arrival time to a congested road preferably in a short form such as elapsed time within a 
mapping cycle, indication on out of mapping range, time related to passing a position such as 
mapping focus, expected time of green light turn on when a road controlled by traffic light 
based on predetermined estimate for the delay in response of vehicle to departure according 
15 to its position in a waiting line preferably in a short form such as elapsed time within a cycle 
such as cycle of mapping samples or cycle of traffic light control (several of such different 
reports can be averaged by the mapping system); reports will preferably use a method of 
transmission that reports characteristic values by transmitting a signal in slot that best 
represents its characteristic value. 
20 Characteristic value of position : 

A value that a probe provides according to a predetermined protocol regarding its 
position, or an indication on its position, such as its distance from a mapping focus along a 
road or otherwise along a path determined by the protocol. 
Ma pping system : 

25 A system comprising a receiver that receives reports from probes and a computer that 

constructs mapping samples from received reports and processes the mapping samples to 
provide characteristics of the road congestion including, but not limited to, one or more of the 
following reports: departure length from the road congestion between mapping samples and 
preferably its varying characteristics; arrival length to the congested road between mapping 

30 samples and preferably its varying characteristics; estimated length of road congestion; 
estimated average waiting time in a congested road; trend in the length of the congested road; 
estimated length of a congested road at a certain time and possibly interpreted values of length 
in a congested road to number of vehicles based on expected average occupation length of a 
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vehicle such as in a stoppage. 

The system will preferably be equipped also with a transmitter that according to a 
predetermined protocol controls the transmission of probes preferably including, but not 
limited to, one or more of the following: required criteria of traffic conditions that enables a 
5 report; resolution of reports; and preferably the time of the transmission of a characteristic 
valueofaposition that relates to earlier time than the transmitted time, disabUng transmission 
of probes that are closer than a certain position to the mapping focus and re-enabling the 
transmission. 

The system will preferably allocate slots to the probes that according to a 
10 predetermined protocol slots divide a range of positions or time interval to smaller segments 
so that each range will be represented by a different slot. 
Map ping sample : 

One or more time correlated characteristic values of position usually relates to time 
constraints that provide a snapshot of probe positions in a congested road. 

15 Kan pe characteristic value : 

A value that represents one or more characteristic values such as positions within a 
range of reports in a mapping sample. Range characteristic values can provide an average of 
positions or average distance from the mapping focus or a weighting average that consider 
parameters that affect inaccuracy in reports. Range characteristic value can also average 

20 several reported values about a common estimate made by probes, for example, estimate of 

green light time setting reported by more than one probe in a waiting line according to 

distance form the traffic light and reaction time to the traffic light. Such reports can use 

differential updates referred to a common time reference. 

Occupation length of vehicle : 
25 ~~ Average segment along a road equivalent to the length between front of one vehicle 

in front or behind of it. 
Mapping ranee: 

A range respective with the mapped part of the road usually covers the congestion 

starting from the mapping focus. 
30 Using the above terminology, in according with preferred embodiments of the 

mvention.mereismusprovidedamemodofestimatmgmeleng^ 

on probe vehicles reporting characteristic values of their position to a receiver of a mapping 
system which processes the reports, the method including: 
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(a) constructing a predetermined number of mapping samples, 

(b) determining in each mapping sample a position that relates to the position of 
a probe relatively far from a mapping focus, preferably a position close to the farthest probe 
position; and 

5 ( C ) selecting from the positions determined in step (b) the position which is the 

farthest from the mapping focus, thereby to determine an indication of the length of the road 
congestion. 

In a preferred embodiment, the position determined in step (b) is the position of the 

farthest probe from the mapping focus. 
10 According to an embodiment of the invention, after construction of a mapping sample 

a response is transmitted to the reporters that disables transmitters that did not transmit a 
report within a preselected range in the constructed mapping sample, to prevent the disabled 
transmitters from continuing to report. The selected range may include the position of the 
farthest probe. 

15 Additionally, in accordance with preferred embodiments of the invention, there is provided 
amethod of determining traffic motion and length of road congestion, in a system wherein 
probes.mresponsetoapredetenninedprotocol, report characteristic values of their position 
to a receiver of the mapping system which processes the reports, the method including: 
(a) constructing a mapping sample that includes at least one of said reports, 
20 (b) selecting a range of said position characteristic values in which the farthest 

reporter from a mapping focus is identified in a mapping sample constructed in (a), 

(c) transmitting to reporters a response that according to a predetermined 
procedure disables transmitters that did not transmit a report within the selected range from 
continuing to report, 

25 (d) receiving further reports and constructing a subsequent mapping sample, 

(e) repeating steps (a) to (d) according to a predetermined procedure, 

(f) selecting from the ranges selected in step (b) the farthest selected range to be 
indicative of the length of the road congestion; and 

(g) determining motion length toward a mapping focus by calculating a range 
30 characteristic value for arange in amapping sample, subsequent to the firstmapping sample, 

which includes the position characteristic value indicative of the closest position to the 
mapping focus and calculating the difference between the said range characteristic value and 
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the range characteristic value of a corresponding selected range in an earlier mapping sample. 

In an embodiment of the present invention, the range selected in step (b) is 
substantially the characteristic value of position of the farthest probe. 

In an embodiment of the present invention, the indication of the length of the road 
5 congestion is determined substantially based on the farthest selected position in the 
constructed mapping samples. 

In an embodiment of the present invention, the resolution of the position reports 
acquired from the probes is determined according to an occupation length of vehicle in a 
congested road in different traffic conditions. 
10 hi an embodiment of the present invention, at least two different reports from at least 

one probe are required to determine length of road congestion. 

In an embodiment of the present invention, the number of mapping samples is 3 to 6 
for expected average percentage of probes in the range of 3 to 5 percent and wherein the time 
interval In an embodiment of the present invention, the number of mapping samples 
15 is determined based on pre-stored data relating to an average motion over time period, 
estimated probability of probe arrival, and estimated arrival rate of vehicles. In some 
embodiments of the invention, the mapping system estimates the probability of probe arrival 
according to a predetermined procedure based on the percentage of probes among vehicles 
arriving at the congestion, and the method may further include: 
20 concatenating aplurality of non-overlapping segments of consecutive mapping 

samples according to motion between mapping samples, 

determining the number of vehicles in the concatenated segments by the ratio 
of the length of the concatenated segments to expected occupation length of vehicles in the 
road congestion, and 

25 determining by a statistical estimator the percentage of probes based on the 

distribution of the accumulated probes identified over the period relevant to mapping samples 
of concatenated segments. 

In an embodiment of the present invention, the statistical estimator is chosen according 
to assumed Binomial distribution of probes in the concatenated mapping samples. 
30 In an embodiment of the present invention, the number of concatenated mapping 

samples is substantially limited to elapsed time interval wherein Ihe expected probability of 
probe arrival to the road congestion is stationary. 

In an embodiment of the present invention, the pre-stored data is optimized to provide 
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thenumberofmapping samples to produce substantially the nunhnuxn expected error in the 
determined indication ofthelengthofthe congestion compared to the real average length of 
the congestion according to the respective mapping samples. 

In an embodiment of the present invention, the optimization criterion is the minimum 
5 difference between the cumulative distribution function of errors that indicates that the 
estimatesaretoolongand^ 
are too short. 

In an embodiment of the present invention, the pre-stored data is prepared based on 
a simulation that shows the number of mapping samples required to provide mimmum 
10 expected error for various conditions of congestion including motion rate, arrival rate and 

percentage of probes. 

In an embodiment of the present invention, at a time prior to determining the 
indication on the length of the road congestion, mapping samples are adjusted according to 
a predetermined procedure that virtually adjusts the position of the mapping focus in the 
15 mappmgsamplesmordertoremovemfferencesbetweenmouonmteand 

In an embodiment of the present invention, at a time after determining indication of 
the length of the road congestion, the determined indication is adjusted by a value which is 
inmcativeofmeprior adjust 
the effect of prior adjustments on the mapping samples. 

in an embodiment of the present invention, successive new indications of the length 
ofmeroadcongestionaredetennmedacoordm^ 

mapping samples. 

Inanembodimentofthepresen^^^ 
to the successive indications on length of road congestion. 

In an embodiment of the present invention, time correlated mapping samples are 
collectedaccordingtoareqmredresolutionofmeroadco 

on departure rate of vehicles from the road congestion. 

In an embodiment of the present invention, the number of vehicles in a lane segment 
of the road congestion is determined according to estimated occupation length of vehicles. 

Inanembodimentofthepresentinvention, the number ofvehicles in a lane segment 
of road congestion is determined according to the estimated occupation length ofvehicles. 

Further, in accordance with a preferred embodiment of the present invention, there is 
provided a method of creating conditions which enable assessment of traffic motion rate 
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toward a mapping focus in a congested road and of the road congestion length at a certain 
time, wherein according to a predetermined protocol probes report characteristic values of 
their position to a receiver of a mapping system which processes the reports, the method 
including: 

5 (a) constructing a first mapping sample that includes at least one of said reports, 

(b) determining a range of said position characteristic values in which at least one 
of said reports was identified in the first mapping sample, and 

(c) transmitting to reporters a response that according to a predetermined 
procedure disables transmitters that did not transmit a report within the selected range of the 

1 0 first mapping sample from continuing to report. 

In an embodiment of the present invention, the characteristic value of aposition is an 
indication of the distance of a reporter from the mapping focus. 

In an embodiment of the present invention, the mapping sample constructed by reports 
is transmitted to the mapping system within a response time synchronized to the mapping 

1 5 system and to the reporters. 

In an embodiment of the present invention, a report of a position characteristic value 
is a signal transmitted by at least one probe in a slot of the response time that is indicative on 
a range of positions. 

In an embodiment of the present invention, the time of the position related reports is 
20 determined by a broadcast query to the probes. 

In an embodiment of the present invention, the selected range in which reporters are 
not disabled includes the position of the farthest reporter from the mapping focus. 

In an embodiment of the present invention, the transmitted response to disable 
transmitters is a message including mapping sample that according to predetermined 
25 procedure reporters disable their transmitters from continuing to report if they had not 
transmitted a report within a range in the mapping sample selected according to the 

predetermined procedure. 

In an embodiment of the present invention, after disabling a transmitter, the probe 
enables its transmitter by predetermined procedure atatime after itpasses the mapping focus. 
30 In an embodiment of the present invention, a non-disabled reporter reports a time 

indication representing its arrival at the reported position. 

In an embodiment of the present invention, the time indication report is a signal 
transmitted by reporter in a slot that is indicative of atime interval in the mapping sample time 
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constraints. 

m an embodiment of the present invention, after a disabling response the mapping 
system receives the new reports and constructs a second mapping sample including new 
reports. 

5 In an embodiment of the present invention, a reporter from the selected range reports 

an indication that it is out of mapping range if it passed the mapping focus at a time prior to 
time constraints of the second mapping sample. Preferably, in an embodiment of the present 
invention, the indication of out of mapping range determines no change in motion rate. 

In an embodiment of the present invention, the indication of out of mapping range is 
10 a signal transmitted in a slot that a transmission in it is indicative on such condition. 

In an embodiment of the present invention, the mapping system determines lengthof 
motion toward mapping focus by calculating a range characteristic value for a range in a 
mappmgsam P lesubsequenttoamsablmgresponse,wmchincludesmeposto 
value indicative of the closest position to the mapping focus and calculating the difference 
15 betweenmerangecharacteristicvalueandmerangecharact^ 

selected range in an earlier mapping sample. 

In an embodiment of the present invention, the mapping system determines the 
departure rate from the mapping focus in a congested road in units of length of a congested 
roadsegmentperumtoftimeaccordingtomedeterminedlengmofmotiontow^ 

20 focus. 

In an embodiment of the present invention, according to motion towards mapping 
focus.memappmgsystemdete^^ 
the road congestion. 

In an embodiment of the present invention, the determined motion length towards the 
25 mapping focus related to average occupation length of a vehicle, including its unoccupied 
space, determines departure rate of vehicle from the congested road. 

In an embodiment of the present invention, the time at which the traffic light system 
changes to a green setting is substantially identified by the mapping system through reports 
from probes. 

30 In an embodiment of the present invention, according to predetermined procedure, a 

probe reports to the mapping system a substantial time when the traffic light system changes 
to the green setting by a predetermined delay taken for a vehicle to react according to its 
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position in a waiting line. 

In an embodiment of the present invention, the time of light change is substantuuly 
identified based on signal transmitted in a slot representing a time interval that best 
characterizes the time report. 
5 fcsomepreferredembodm^ 

procedure, the mapping system estimates arrival rate to the congestion, and the method 

further includes: 

concatenatmgapluralityofnon-overlappmgsegmentsofconsecutivemappm^ 

samples according to the two said position related reports, 
! o determining time intervals between two time of arrival reports corresponding 

to two position related reports, 

determining average arrival rate in terms of segment length occupied by 
arriving vehicles between successive mapping samples by calculating the ratio of the total 

lengmof concatenated seg^^^^ 
15 In an embodiment of the present invention, the number of concatenated mapping 

samples is substantially limited to elapsed time interval in which the average arrival rate of 

probes to the road congestion is expected to be stationary. 

In an embodiment of the present invention, the length of motion of non-disabled 

reporter towards mapping focus between two consecutive mapping samples determines the 
20 point between consecutive concatenated segments. 

RPTF.y TreSCRIPTTnN r>F THE DRAWINGS 
In ordertobetterunderstandtheinvention and to see how the same may be carried out 
in practice, non-limiting preferred embodiments of the invention will now be described with 
reference to the accompanying drawings, in which: 
25 Fig. 1 shows an initial map generated in a prior art ITS system; 

Fig. 2 shows a second, more detailed map, generated during a second iteration in the 

prior art ITS system; 

Fig. 3 shows a graph of additional information which is generated in the prior art ITS 
system; 

30 Fig. 4 shows a graph of further additional information which is generated in a prior art 

ITS system; . 
Fig. 5 is a general block diagram of a transmitter for the prior art ITS system, which 

SUBSTITUTE SHEET (RULE 26) 



PCT/IB00/00239 

WO 00/54143 _|g_ 

is also useful in the present invention; 

Fig. 6 is a block diagram of a receiver for prior art ITS system, which is also useful 

in the present invention; 

Figs. 7-12 illustrate a scheme for traffic jam mapping, in accordance with a preferred 

5 embodiment of the invention; 

Fig. 13 shows a system for acquiring mapping information, in accordance with a 

preferred embodiment of the invention; 

Figs. 14 and 15 are schematic block diagrams of traffic jam mapping systems in 
accordance with preferred embodiments of the invention, as interfaced with a car navigation 
10 system; 

Fig. 16A shows exemplary results of computer simulation of queries of vehicle 
positions with respect to a focus of congestion, in accordance with a preferred embodiment 

of the invention; and 

Fig. 16B shows a stream of traffic constructed from the queries of Fig. 16A; 
15 Figs. 17Aand 17B are a flow chart of a method for determining the optimum number 

of maps for estimating the length of a queue; and 

Figs. 18A and 18B are tables showing computer simulation results of application of 

the method illustrated in Figs. 17A and 17B. 

TYFTATT ED PESO* TPTTON OF A PRFFFPFFP EMBODIMENT 
20 In large scale road traffic mapping systems, it is generally desirable to have a large 

enough pool of transmitting vehicles in the region such that there be a good chance that at 
least one and preferably a plurality of vehicles with transceivers transmit, according to 
mapping queries. The percentage of vehicles that are equipped with such transceivers 
determines levels of probability of success of identification of a stoppage or slowdown in 
25 traffic and the level of mapping the congestion. 

In the embodiments of vehicle traffic mapping systems, described in conjunction with 
Figs. 1-4, a focus of the mapped congestion is identified without any need for a priory 
knowledge of intersections (which may be a focus of congestion) by the broadcast query 
system. The efficiency of this strategy can be improved, in a preferred embodiment of the 
30 invention, by using maps within the broadcast system in order to identify and determine a 
tentative focus of a congestion, as for example an intersection. In addition, no mater how a 
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focus is determined, the efficiency of the system can be improved by vehicles reporting their 
position (inresponsetoaquery)as distance and direction fromafocus, or ifthey have internal 
maps, as distances along a branch from an intersection or focus of a disturbance on an open 
road. This can provide a relatively high accuracy of estimated length of jams based on 
5 distances from the focus, utilizing a relatively smaller percentage of vehicles being equipped 
with mapping transceivers. Preferably, this will utilize and assume a computation logic that 
logically fills the gaps between the focus and the measured length by artificial responses of 
vehicles on the (virtual) map, experiencing the same congested conditions. Effective use of 
repeated mappings, as described below, can vastly improve both the probability of detection 
1 0 of congestion and the accuracy of its mapping. 

Alternatively, the focus can be determined at a central station and broadcast to the 
vehicles for use as the reference position for reporting distances. 

In a large city, reliable mapping might require a large number of vehicles, albeit a 
small percentage of the total number of vehicles equipped with mapping transceivers. 
15 Furmermore,evenwhenoneofmevehiclesismaqueueitismfficulttoestimatemep^ 

of the end of the queue from an instantaneous position. 

In accordance with a preferred embodiment of the present invention, jam length 
mapping of a slowdown or stoppage is based on a distance and direction of vehicles from a 
focus of a slowdown. This distance may be an intersection or, in the case of an accident, for 
20 example, it may be the site of the accident. Use of distance and direction (initially) from a 
focus as foe fransnuttedinfo^^ 

for position of the vehicle, is more bandwidth efficient since transmission slots need be 
provided only for one dimension of positions rather than for a two dimensions. 

In a preferred embodiment of the invention, where a vehicle reports a slowdown or 
25 stoppage (as a result of a query such as described in conjunction with Fig. 2 or 3, or some 
other method), while traveling towards an intersection, the focus is assumed to be at the 
mtersection.Formesefoci,mequeriesandpositionsthatarebroadcas^ 
from an intersection. If no intersection appears to be related to a slow down, the focus is 
unknown, possibly until a first vehicle passes the focus and speeds up and can respond to a 
30 query that maps such a condition. In order to identify the focus and confirm the focus as an 
intersection when it has been assumed to be such, a vehicle preferably records the time and 
position at which its speed increases to indicate that it has passed the focus. Then query based 
or other methods maybe used to identify andmap one or aplurality of foci. For example, such 
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a query could include the criteria that vehicles which were in the slowdown and speed up (at 
afocu^respondmoneoraplu^^ 

positionontheroadatwhichthey speedup and optionaUy, time indicative slots corresponding 

to the time and position at which they speed up. 
5 suchqueriescanbeextendedtogiveanindicationofmepositionofmefocusandthe 

length of time it takes a vehicle to traverse the slowdown. This can be implemented by the 
same or another query that use the current timeasatime reference in the query, with the time 
at which the vehicle first experienced the congestion being the requested information in 
respectiveslots.Thus,mefocusandelapsedtimeexperiencedby^ 
10 provides an indication of the traverse time for the congestion. 

In a preferred embodiment of the invention, means are provided for minimizing the 
false alarms in slots that cover the portions of roads beyond the actual congestion. One 
strategy to minimize such false alarms is to minimize the length of the road, leading to the 
congestion, which is mapped. Under this scheme, the number of slots (i.e., the length of the 
15 roadbeingmapped) is gradual 

term statistics of the current mapping process. Thus, for example, the length could be some 
multiple of the distance from the focus of a previously identified vehicle, pass-out speed of 
the vehicles, etc. Another way of reducing false alarms is to correlating successive mapping 
cycles to determine "non-repeats." Thus, for example, if on one mapping cycle a vehicle at a 
20 given distance from the focus reports a slowdown, this indication would be ignored if on a 
successivere-mappmgc^ 
the vehicle could have progressed, is reported. 

Alternatively or additionally, the position of the focus of the slowdown is identified 
fromasource external to the systemorby any othermethod. Before the position ofthe focus 
25 is known it may be estimated as the position ofthe most forward vehicle in the congestion. 

In a preferred embodiment of the invention, length measurement from a focus (in 
which a focus has been assumed or positively identified) is used for periodic remapping. In 
such a system, repeated remapping queries are periodically broadcast requesting that slowed 
vehiclesmovingtowardasame focus to transimtmek positions wim respect to distance from 
30 the focus. This will continue as long as there is an indication of congestion around the focus. 
The data from successive cycles is used to fill in the interpreted picture of the remapped 
congested road, based on slow changes in the position of vehicles and vehicles that pass into 
and out ofthe congestion. This may be achieved by mapping the congestion periodically and 
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combinmganumberofsuc^ 
somewh^lower.butstiU^ 

broadcast in a number of slots in the same mapping cycle, where they broadcast in all slots 
which represent their position within a given prior period. In both cases, each vehicle 
5 represents (virtually, for the map) a number of vehicles at different positions withm the 

slowdown. 

Also,forfocusesmataremte^ 
for example when using query based periodic remapping, to improve the accuracy of the 
remappmgbysynchromzingitwimmeperiodicclearmgsofmeconge^ 
L0 caused by traffic light changes. The synchronized time could be identified by focus/tune 

related reporting as described above. 

In preferred embodiments of the invention, it is not necessary for the portion of the 
system that is in the vehicles to have reference maps. Rather, the queries request information 
from vehicles traveling toward the focus, that are in positions of interest. 
15 Altemativelyoradditionally.mevemclemayassociateitspositionwim 

inaroadmapmemory.Underthese^^ 

thepurposesofinitialqueriesto local groups and sub-groups, corresponding to areas. Dunng 
mapping, when the resolution reaches a given level, or during remapping, short map related 
codes based on road coordinates and on the position of the focus, are transmitted accordmg 

20 to assigned "one dimensional" slot allocation mapping. 

Anindicationfromasmglevemclematitisstoppe4orslowedatsomedistancefrom 

anmtersection(orfromsomeo^ 

to an actual traffic stoppage. For example, the car may be stopped or slowed due to a non- 
tmfficcause such as a malfimction which may neither affecttrafficorbetheresultofatraffic 
25 condition. Furthermore, if the proportion of cars in the total population is small, the 
probabiUtytluitacarwillbeatornearmeendofmetrafficslowdownma^ 
high Thus, the measurements of the length of the stoppage or slowdown have an mherent 
accuracy depending on the percentage of vehicles with mapping transceivers in a region and 
in a mapped road. It is desirable to increase the resolution of the measurement of the length 
30 of the stoppage or slowdown and to increase the reliability of determination that a stoppage 
or slowdown exists, while using a relatively small percentage of transmitting vehicles. 

In general one or more of the following means and methods for improving the 
rehability and/or r^ 
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that are equipped with intelligent mapping transceivers, are provided, in accordance with 
preferred embodiments of the invention: 

1) Requiring that at least two distinguishable vehicles be stopped at a same traffic 
stoppage, otherwise ignoring the transmission that represents a single vehicle (i.e., a single 

5 slot). 

2) Where there is slow moving traffic (with vehicles passing through and out of the 
congestion and new vehicles passing into it), estimating the end of the traffic slowdown by 
the furthest vehicle from a focus of a slowdown (an intersection or an accident, for example). 

3) Repeating the mapping during aperiod of time that is longer than a relatively short 
10 mapping cycle time and providing a map based on two or more of such repeated cycles. This 

multiple remapping information may be transmitted during multiple cycles or during a single 
cycle as described above. Since one can expect (and predict statistically) that additional 
transmitting vehicles will enter the congested road, this remapping method gives a better 
estimate of the length of a slowdown, albeit with a loss in time resolution. Note that this may 
15 result in the same vehicle being mapped during multiple remapping cycles or it may result in 
different vehicles being mapped from different cycles. This method is based on the 
assumption, which is usually correct, that within a one or a few minutes using two or more 
mapping cycles there is no significant change in the length of the queue on the congested road. 
Thus, the loss of time resolution of using multiple remappings in mapping a jam is not 
20 significant. It is possible that, if too long a time, or too many repeats are used in forming a 
map, the map may represent a statistical peak in the length of the congestion. This can be 
avoided, if necessary, by limiting the number of multiple queries utilized in forming a map. 

In preferred embodiments of the invention, there is an overlap between successive 
groupings of mapping cycles (i.e., windows) that are used to form an estimate of the length 
25 of the congestion. That is to say, in a subsequent estimate, the oldest mapping cycle result in 
the previous estimate is discarded and replaced by a later mapping cycle result. In this way, 
the result is updated more often than it would be if non-overlapping groupings were used for 
sequential estimates. Alternatively, non-overlapping groupings are used. This results in 
correlation between successive windows. However, to avoid occasional large sharp 
30 fluctuations in the queue length, it may be desirable to filter the output stream with a filter 

such as a median filter. 

4) Analyzing the time development of the motion of a vehicle farthest from the 
intersection or other focus of a slowdown. Since vehicles periodically enter the queue, a plot 
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of the position of the farthest vehicle will have an undulating, if irregular, form, with the 
greatest distance (peak of the undulation) representing the position of the end of the queue. 

5) Evaluating the reliability of responses from vehicles substantially farther from the 
focus than a previously determined end of the congestion based on pass-out rates and on an 
5 expectationofanadditionalresponseinasimilarposition(orsomewhatcloserposition,based 

on an estimated vehicle velocity through the congestion) from the vehicle. 

In some preferred embodiments of the invention, each vehicle keeps track of its last 
few positions and transmits, as its position, its largest distance from the focus of the 
congestion, while it is moving slowly. In other embodiments, each velucle transmits its actual 
10 relative position (preferably in a slot in response to a query) and either a central station or 
other vehicles, as described below, calculate an estimated length of congestion. 

For clarity of presentation, the following definitions of terms are provided: 
l)a"map"isasingleresponsetoafocus related mapping query, generally based on responses 

in allocated slots; 

15 2) an "overlay" is one of a series of maps that is used in detenrnning an estimate of a length 
of queue in a congestion situation; 

3) overlaying means estimating the length of a queue by overlaying a series of maps such that 
the longest queue estimate from these maps (i.e., the response that is farthest from the focus) 
is used as the estimate of the average length of the queue; 
20 4) a "window" is a period during which overlays are combined to estimate an average length 

of a queue; 

5) a "probe" is a vehicle that can respond to system queries; 

6) a "mapping cycle" is the time between two overlaid maps; 

7) "departures" is the number of vehicles that clear the queue, past the focus of congestion, 

25 in a given cycles; 

8) "average departure rate" is the average number of vehicles that clear the queue, past the 

focus of congestion, per cycle; and 

9) "average arrival rate" is the average of a Poisson Probability distribution that governs the 
arrival of vehicles at the end of the queue. 

30 In general, it is applicable to estimate the average length of the queue, when the 

statistics governing the process is stationary. When the average queue length is changing (as 
in developing congestion) it may be desirable to report either the average length of queue 
during the multiple queries or the (statistically) average queue length referenced to the last 

SUBSTITUTE SHEET (RULE 26) 



PCT/IBOO/00239 

WO 00/54143 24- 

overlay in the window. The term "average" queue length as used herein is the queue length 
averaged over a number of trials (e.g., a number of overlays), to reduce the effects of 
statistical variations in the length of the queue. Such statistical variations occur in any 
statisucalprocess,evenifmeprocessisstationary,i.e.,itdoesnothave^ 
5 if the process is not stationary, i.e., it does have a trend or other determined variation), it is 
generally desirable to first remove the effects of the trend or determined variations, for 
exampleusingoneoffcemethods^^ 
to estimate the center of the variations. 

As shown below, the main statistical process is the number of vehicles entering the 
10 queue during a mapping cycle, which has Poisson statistics. When the clearance rate of 
vehicles per cycle is a constant, the length of the queue also has Poisson statistics, which may 
be stationary or not depending on the relationship between the average rate of entry of 
vehicles and the clearance rate. When these rates are the same, the statistics are stationary. 
However, this is seldom the case in the real world, especially on a short term basis and thus, 
15 in order to evaluate the queue length with best accuracy, it is desirable to reduce all non- 
Poisson statistical effects on the queue length. 

It should be noted that just as too few queries (e.g., when the queue is too short) will 

resultmerrcrsinmelengtho^ 

may also cause errors, although for a different reason. As with any statistical process, if the 
20 furthest distance is based on a large number of cycles, and the length within each cycle is 
determined to be the largest distance of all the samples in the cycle, there is an increasing 
statistical probability that the length will be overestimated by an overlaying process. The 
statistical probability for such errors will depend on a number of parameters, namely, the 
percentageofprobesmmetotalnumberofcars,menumberofqueriesmatare»^ 

25 other parameters described below. 

Forasmallpercentageofprobes,thestatisticsarepoorerandmeovemheln^ng 

for small numbers of overlaid maps there is a relatively high probability that no or too few 
cars will be in the queue and that the end of the queue will be poorly determined for this 
reason. In order to assure that the queue be optimally mapped, a number of maps should be 
30 overlaid,uptomepointwheremeprobabilityoferrorduetotoolonganestimateofthequeue 

length is as large as the probability of the queue being measured as being too short. 

On the other hand, when a larger percentage of probes are present, the probability of 
under reporting the queue length is reduced fairly quickly with number of overlaid maps, 
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while the probability of over reporting is relatively high with increasing number of overlays. 

It should be understood that the change in expected statistical errors caused by using 
a larger or smaller number of overlays for a single queue length estimation is generally not 
very large, under many typical conditions, if a reasonable number (such as four or five) 
5 overlaid queries are used. It has been determined by the present inventor that such typical 
conditions may be present if the percentage of probes is between 3-5%. However, for greater 
accuracy and under certain circumstances, the accuracy can be significantly improved by 
optimizing the statistics of the queries (i.e., the number of overlaid queries) and by 
synchronizing the queries with some cyclic event (such as the changing of a traffic light). The 
10 present inventor has performed a number of statistical studies to determine the optimum 
numberofoverlays,utilizingthe first criteria described above. The numberof overlays varies 
between about 3-5 for large percentage of probes (such as for example 3-5%) to 7 or 8 for 
small percentages of probes (e.g., approximately 2 percent), when a road with two lanes is 
joint, for mappingpurposesasasingle double density lane. Not unexpectedly, the mean error 
15 mestimatedqueuelengmislargerforsmallpercentagesofprobeswanforlarg^ 

The parameters that appear to affect the (statistically) optimum number of overlaid 
queries and the expected errors for various situations have been determined by the inventor. 
In making this determination, the inventor has performed statistical studies of the (assumed) 
actual and reported queue lengths. In these studies the following factors have been found to 
20 affect the number of overlays to be used for statistically optimum results and of the 
statistically expected errors: 

(a) the percentage of probes reporting in a sampled map; 

(b) the average arrival rate per map (stationary Poisson distributed); and 

(c) the departure rate per map (or the actual departures if known). 

25 Furthermore, these studies have shown that synchronizing the timing of the queries 

with the traffic light, reduces the effects pf cyclic variations on the estimation of the length 
of congestion and assists in determining the average occupation length of a vehicle in the 
queue, e.g., based on map density, thereby to convert the size of the queue from length units 
into number of vehicles. This may be useful in estimating departure rates, estimating arrival 

30 rates and estimating the percentage of probes, as described below. Since it is not always 
possible to query the vehicles in the queue at exactly the desired time (for example, when two 
congestions, having the same query time are being simultaneously mapped using the same 
communications resources), in some preferred embodiments of the invention, the vehicles 
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store position information, and provide this information to a later query. The vehicles may 
store position information for a specific time (based on a previously given command) or they 
may continuously store such information. 

Thus, in order to optimally map congestion it is useful to know all of these factors in 
5 determining the number of queries to overlay in a mapping cycle. 

mapref OT edembodimentofmemvention,mestatistica% 
to be overlaid is determined based on the statistical probability of different types of errors. 
The criteria for determining the balance between the errors which is considered optimum is 
in some sense arbitrary, and is based, inter alia, on an assessment of which type of error is 

1 0 more problematic. 

One possible basis for a criteria for determining the number of maps to be overlaid is 
to compare the probability that there will be an overestimate of the length to the probability 
that the system underestimate the length of the queue. One possible criteria is to set the 
number of maps in a window equal to that at which the two probabilities (and the average 
15 length ofthe error) are equal. Fig. 1 7 shows a flow chart of a method, utilized by the present 
inventor to determine the statistically optimum number of overlays. Figs. 1 8A and 1 8B show 
a table containing some ofthe results of this study, for average departure rates of 20 vehicles 
per cycleandlOvehiclesper cycle, respectively. It should be noted that the flow chartofFig. 
17 is limited to a typical situation in which the average departure rate is equal to the average 
20 arrival rate. Figs. ISA and 18B show the results of a more general simulation in which the 
average departure rate is not necessarily equal to the average arrival rate. 

Other methods of determining the optimal number of overlays may be used. For 
eX ample,insomesituat^ 

noted, In such situations, the number of overlays may be increased. In others, fast response, 
25 rather than accuracy may be important which would lead to a reduction in the number of 
overlays used. 

As indicated above, the determined optimum number of overlays is dependent on a 
number of parameters, one of which is the percentage of probes in the total pool of vehicles 
and, more importantly at, or in, the queue whose length is being measured. While one can 
30 knowmeproportionofprobeveWdesmmetotalveMclep^mation.itU 

a priori, the proportion of probe vehicles on the road as compared to the total number on the 
road, at any one time. It is even more difficult to estimate, a priori, the proportions at any one 
site, at any one time. It has been found in statistical studies carried out by the inventor, that 
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having 3%-5% of probes is generally sufficient for giving fairly accurate estimates of queue 
length. Lower percentages can also give useful results while larger percentages give more 
statistically accurate estimations of queue length. Providing such a large percentage of probes 
in the general "population" of vehicles can be very expensive. However, one possible way to 

5 provide for a wide distribution of probes (without the great expense of equipping large 
numbers of vehicles that are used only sporadically) is to equip vehicles that are on the road 
a large percentage of the time, such as taxis and/or buses and/or other types of commercial 
vehicles, with transmitters. However, the distribution of such vehicles is not uniform (poor 
neighborhoods have few taxis present, for example) and varies with the time of day and the 

10 influx of vehicles during rush hours (during which times automobiles from outside the area 
increase the total number of vehicles without changing the number of probes). 

Thus, while models may be constructed for detemining an optimum number of 
overlays based on any particular criteria and situation, the required knowledge of the actual 

parameters is unavailable. 

15 In a preferred embodiments of the invention one or more of the parameters affecting 

the statistically optimum number of overlays is estimated based on either "hard information" 
(direct information with respect to the changing times of traffic lights, for example), "soft 
information" (such as statistics of probe percentages based on prior surveys) or by the 
estimation of these parameters from the responses to the queries, themselves. 

20 Fig. 16Aisarepresentativeexampleofaseriesofresponsesto successive queries that 

aremadeatasourceofcongestion. It is present^ 

preferred methods of estimating the various parameters necessary for determining a 
statistically optimum number of queries (overlays) to be used in a mapping. 

As a first step in this estimation process, the number of vehicles passing through the 

25 congestion (vehicles per mapping cycle) may estimated using any suitable method. Aside 
from its importance in detenniiiing the optimum number of overlays to be used in estimating 
the length of the queue, accurate determination of velocity can also be an important 
consideration in traffic control applications. In particular, the measurement of velocity is 
important in determining the delay time at a particular intersection and the average number 

30 of vehicles passing the intersection per/mapping cycle. The average delay time will be an 
important consideration in decisions related to rerouting of vehicles and the average number 
passing the intersection (in relation to the length of the queue) is useful in determining an 
optimum division of the cycle time of the traffic light among the different streams of traffic 
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entering the intersections. However, due to the nature of movement of vehicles in congested 
situations, accurate measurement of velocity should preferably be made with care. 

The motion rate of vehicles through the congestion can be estimated in several ways, 
in accordance with preferred embodiments of the invention. 
5 A first method of determining the motion rate of vehicles is to compare patterns in 

successive maps. Since successive maps show substantially the same patterns, displaced 
toward the focus of the congestion by the motion of the traffic, the movement per mapping 
cyclemay be determined simply bedetermined from the movement of the pattern toward the 
focusofthecongestion.Suchpattems 
10 in adjoining lanes, since such different velocities will result in small changes in the patterns. 

In accordance wim an alternate preferred embodiment ofthe invention, the movement 
per mapping cycle is determined by determining the distance traveled by a selected vehicle 
between queries. Since the vehicles bear no identification, in preferred embodiments of the 
invention, this requires singling them out, based on some characteristic of the vehicles. One 
15 way of doing this is to find a relatively isolated reporting vehicle in the queue, estimate a 
velocity and look for the vehicle at the next query at a position estimated from the average 
velocity. This may be done with vehicles in the first few clearance lengths from the 
intersection or far from the intersection. A clearance length may be defined as the length of 
the queue that clears a traffic light for each cycle. This method may be advantageous because 
20 it may not require any special queries or calculations by the vehicles themselves. 

Alternatively, the vehicles may be asked specifically to respond or not to respond, or 
to determine their own average velocity or distance traveled per cycle and asked to transmit 
this value during designated slots in a query. 

Alternatively, in a special query, only the last vehicle in the previous map may be 

25 askedtotransmititspositio*^ 

that the distance traveled may be estimated. 

Alternatively, a special query may ask only a particular vehicle (based on its previous 
position)tobroadcast during thequery.This enables isolation ofthe vehicle even whenmany 
probes are present. It is often simple to request that only the last probe vehicle in the queue 

30 broadcast in the following map. This vehicle will then become the reporting vehicle closest 
to the focus in the next map. In order to minimize the number of special queries, a 
predetermined protocolmaybe used whereby only the farthest vehicle from the focus may be 

iratmctedtorespondina^^ 
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may broadcast their positions as usual. 

These last WO methods are illustrated in Fig. 16A, wherein vehicles are shown as 
boxes inaqueue, each box representing a slot related toamapping length of a vehicle,™* 
movement being to the left. Empty boxes represent a vehicle space, occnpied by an average 
5 hmgth non-reporting vehicle and boxes with numbers represent reporting vehicles. In both 
simulations and actual mapping queries, as described herein, distences are denomtnated m 
"average vehicle lengths" or any other predetermined occupation length of a vetacle m a 
queue. It should be understood that reporting is based on distance from .he focus, which for 
computational purposes may be replaced by an average vehicle occupation length. The 
10 average occupation length may be defined as the average occupation length it, a standtng 

take into account a varying queue density. 

Since Fig 16A represents a simulation, the actual nnmbers of vehicles in the queue 
a, any one time is known (to the computer performing the simulation. Thus the "reported" 
15 ^ofmesmdyc^mm^comparedtou.e-actiud-valnes.I.shonldbeno.edhowever, 

^ ^.«*»^^*w-^»^*^ ,- * ,,l,,Hh,-w ■■ , * 

,6A shows a fairly constant spacing and throughput rate. This is usefitl for exposrtory 

purposes, but is not generally the case. 
2 „ From Fig. 16A, it is clear that vehicle 17 advances 10 predetermined velucle 

between tire second rctd third queries; and vehicle 37 advances 10 lengths between the nurd 
,„dfburth,»erie*.Howev.r,,woproblemspr«ent thrives. Fir*, in me actud sftuaton, 
fte vehicles are no. identified. Second, no estimate of movement can be made between maps 



25 4 and 5. 



Asm&atedabove.ooewaymovercome.hefire^ 

of the invention, is to perform a social trial (query and map) in which only the iartites, 
vehiclemmequeueformepreviousqueryfandpo^iblynewvemcleswhichjoinntequeue), 

determined according to the responding slot which is known to the responded is asked * 
30 ^nd.totHsway.thefira.tc.oses.^ 

identified as me last vehicle in the previous map. Thus, its movement during tire mappmg 
eyele (and the movement of tit. queue as a whole) can be estimated. If new vehicles also 
respond to the query to form the new map, there is no loss of information needed for the 



SUBSTITUTE SHEET (RULE 26) 



„ mMtm PCT7IBOO/00239 
WO 00/54143 3() ^ 

determination of queue length by performing this type of query. 

With respect to the second problem, in situations where the farthest probe moves out 
of the queue, no measurement of throughput is made. Unless the percentage of probes is very 
small and the queue length is very short, this should not usually be a problem, since if it 
5 occurs the estimated average motion per cycle can be used to estimate and bridge the gap. 

Inapreferredembodimentofthe invention, the percentage of probes is estimated from 
the actual data, utilizing a concatenation of data from successive maps. 

In order to demonstrate how this is done, a single stream of traffic is constructed from 
the data of the five maps is shown Fig. 16A. This single stream is shown in Fig. 16B. As 
10 indicated, since the position of the reporting probes in the stream of traffic is known, these 
positions can be used to act as a bridge connecting the maps whenever a probe is present in 
two succeeding maps. This allows a simple connection between maps 1-4 

With respect to connecting maps 4 and 5, since an average throughput is known from 
thepreviouslydescribed estimation, the position ofthe first vehicle ofmap 5 can be estimated 

15 asbcangthevdricleatoto 

sole vehicle to report in map 5 is placed fourteen slots after the vehicle reporting in map 4. A 
similar estimate may be made even if no vehicles respond to a given intermediate query. 

Based on the density of probes in this stream of traffic, the probe density is estimated. 
The distribution of distances between vehicles is a geometric distribution Since the distance 
20 betweentwoprobesismemeasuredlengmdi^^ 

then the number of vehicles between probes may be easily estimated. The probability of 
probes can then be estimated (for example using maximum likelihood methods), as 1 per 
mean distance between probes (in vehicles units). To improve these statistics, the length of 
the concatenation should be as long as possible and may include maps which span more than 
25 a single window and which include may past maps. 

A third parameter that is known to effect the number of queries to be used in a 
mapping is the estimated arrival rate of vehicles. This can also be estimated from the stream 

of traffic shown in Fig. 1 6B. 

In order to estimate the arrival rate of vehicles per mapping cycle, time related 
30 positions on the stream of traffic of Fig. 16B should first be determined. Then the number of 
arrivals per generated stream oftraffic is estimated. Finally, me arrivd rate per mapping cycle 

can be estimated. 

In order to determine the time of arrival (vehicles/time) a first vehicle and its arrival 



SUBSTITUTE SHEET (RULE 26) 



PCT/IB00/00239 

WO 00/54143 _ 31 _ 

at the far end of the congestion is determined. Then the arrival time of a second vehicle at the 
far end of the congestion is determined. The number of vehicles that arrive between the time 
of arrival of the first and second vehicles is estimated, based onthe stream of traffic (Fig. 16B 
andtheaverageveMcleoccupat^^ 
5 from the ratio between the number of vehicles (occupation lengths) that arrive between the 
first and second vehicles divided by the difference between the times of arrival of the two 
vehicles at the congestion. The average number of arrivals per cycle is then determined by 
multiplying the average arrival rate by the map cycle time. 

In a preferred embodiment of the invention, vehicles note their time of arrival at a 
10 congestion. When the arrival rate is to be calculated, a particular vehicle (which can be 
identified from its position in the current queue) is instructed or allowed to report its arrival 
time at the congestion, for example in special slots. A stream of traffic (such as that shown 
in Fig. 16B for the computer simulation) is constructed, for successive queries preferably 
utilizing the same methodology described above: When the stream is long enough, a second 
15 vehicle, identified as the one near the end of the stream in the previous map, is requested to 
transmit its arrival time at the congestion. This vehicle is also identified from its position in 
fcemencurrentqueueasmefiKtrespon^ 

In a preferred embodiment of the invention, the number of vehicle lengths used in the 
queries andmapsadapts to expected ohanges in the length ofthe queue. In one methodology, 
20 the length ofthe queue is based on previous estimated lengths and on any trends noted in the 
length. Alternatively or additionally, the length ofthe query (number of slots) depends on the 
net estimated arrival rate of vehicles, which are expected to result in a change in the length 
ofthe queue. 

In each ofthe above determinations, a single lane of traffic is assumed. If multiple 
25 lanes of traffic are present, the maps may or may not differentiate between the lanes. If they 
do not, then the calculations do not differentiate between the lanes and the stream of traffic 
may be considered to be a "double density" stream. If the lanes are mapped separately, the 
data can be combined and the traffic flow calculated on the basis of double density. 
Alternatively, the congestion in each ofthe lanes can be determined separately. 
30 When multiple lanes are present, special considerations sometimes have to be taken 

into account, especiaUy for certaintypes of congestion. Itparticulax,inprefe^ 
ofthe invention, the effects of turning lanes and the effects of merging traffic in a common 
direction are taken into consideration. 
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In accordance with a preferred embodiment of the invention, where turning lanes or 
merging traffic (as at a repair or accident site) are present, it is normal to expect traffic 
velocities to be different for different lanes in the same direction and timing of traffic s,gnals 
to be different for the different lanes which move in different directions. One solution to the 
5 problemsmusrdsedistoqueryea^^ 

communications resources. Another possible solution is use information for determining the 
motion rate only from vehicles that are sufficiently far from the focus of the congestion such 
that the effects of the turning lane or merge are "homogenized, i.e., for which the effects of 
the congestion are the same for all the lanes.. It is also possible to instruct vehicles which are 

10 turning not to respond to the query or to respond to a separate query. 

Momerprobleminaccuratedetermmationofmelengmofthequeueismatm^ 

undulates depending on the phase of a traffic signal (where the queue is caused by simple 
congestion at a crossing). As indicated elsewhere, it is possible to measure such undulations 
anddetermineavaluebasedonapa^ 
15 number of measurements over a long period of time, since not only must the undulations be 

found,butmeyarepreferabty^ 

traffic lights is known, as for example where the system is part of a larger traffic control 
system, this information can be used to synchronize the measurement with the cycle of the 
traffic lights, which presumably controls the undulations in the length of the queue. Use of 
20 suchsynchronizationallowsformeremovaJofmeundulations 

valuesofthequeuelengm^ 

easily detenmned by mappmg slots representing relevant!^ 

average occupation length of vehicles. 

Alternatively, the timing of the traffic light can be determined automatically from a 
25 special query of vehicles which pass the traffic light. If such vehicles transmit their time of 
passing of the traffic light, this time can be combined with an estimated time for the vehicle 
totravelfromitslastknownposto^^^ 
relationsmpoftime^ 
to green. 

30 There are several problems with using synchronized systems to measure length. One 

problem is that the effect of the traffic light change travels down the queue at some average 
velocity Thus, while vehicles near the traffic light will begin the move soon after the light 
changes to green, andstop soon after the light changes to red, vehicles further away may only 
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start moving only after the light has turned red. In extreme situations, vehicles may move in 
response to a previous cycle of the light and not in response to the current one. In many 
instances, vehicles far from the traffic light may actually move slowly all the time or be 
subject to stops and starts that are unrelated to the cycle. 
5 Since there is much interest in measuring long queues, the measurement of such 

queues must be undertaken with care when using a synchronized system. In a preferred 
embodiment of the invention, the queries are synchronized with the traffic light cycle, 
preferably, with the turn on of the green (or just prior to the turn on of the green). However, 
the response of the vehicles is based on their own experience, such that not all the vehicles 
10 moveinsynchronism.Forexample,arelativelyoptimaltimeformati^^ 

the light changes to green, because at this time the longest line of vehicles will be stationary 
and the predetermined average occupation length of vehicles can be more easily and 
accurately predicted (e.g., for the purpose of estimating the queue length in terms of number 
of vehicles). When such a query is made, vehicles which are stationary report their present 
15 position. Those vehicles that are moving in a cyclic fashion report their position the last time 
theywere stopped. Those that are movingslowly report their actual position at thetimeofthe 
probing. This is believed to give a stable and consistent measurement of queue length. It 
should be understood, that due to instantaneous radio communication band-width limitations, 
the probing time may be different from the transmit time due to delays in sending the query 
20 or the responses, i.e., a query may ask for responses based on conditions of a previous map. 

In accordance with a preferred embodiment of the invention, a series of computer 
mapping trials are performed (for example in accordance with Fig. 17) in which the errors 
with time are computed as a function of the various parameters given above. In a preferred 
embodiment of the invention, when mapping congestion, the parameters necessary for 
25 detemimng me optimal numb^^ 

queue length are determined based on the mapping trials. As set of the optimal number of 
maps are then combined to perform the mapping of the average length in a given time 
window. It should be appreciated that, while the above described methods are preferred for 
calculating/estimating the desired parameters, other suitable methods may also be used for 
30 calculating/estimating these parameters . 

In a preferred embodiment of the invention, the longest distance from the focus to any 
probe for all the maps in a window is determined as being the estimate of the average length 
of the actual congestion. However, if the average rate of clearance of vehicles per cycle from 
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the congestion is substantially different (e.g., more than about 5%) from the average rate of 
arrival of vehicles per cycle at the congestion, the distance of the last vehicle in each map is 
preferably artificially adjusted to account for the trend in the length of the congestion caused 
by the difference in arrival and clearance rates. The object of this adjustment is to remove any 
5 trend in the statistical distribution so that it is substantially a Poisson distribution and, then 
to estimate the average length of the congestion as described above and, finally, to readjust 
the estimate to correct the artificial adjustment 

One wayof adjusting to remove the influence of trends on the average arrival rate and 
irregularities from the departure rate is to use the concatenation of the maps in a window and 
10 to reconstruct maps that do not include irregularities in the departure rate per cycle or trends 
in the average arrival rate. The departure rate per cycle as well as the average arrival rate can 
be calculated using the processes described above. This adjustment is preferably executed on 
the departure side of the road congestion, e.g., to compensate for a generally extending road 
congestion, the departure rate may be increased to shorten the mapping sample at the 
15 departureend. The adjusted maps in me window can then be used, in an overlaying process, 
to estimate the average length of the queue in the window., e.g., by selecting the farthest 
respontingposition(wimreferenceto^^ 

of the window may be selected to include the number of maps that provides a minimal error 
in estimating the desired parameters. The tables in Figs. 1 8A and 1 8B indicate, inter alia, the 
20 optimalnumbersofmapsinamappingwindowofasy S temusing4percentprobes anda 
system using 3 percent probes. It should be noted that the percentage ofprobes (which is one 
of the parameters affecting the window length) may also be estimated, initially, and is 
consequently adjusted based on a concatenation of several maps before determining an 
optimal mapping window. Readjustment of the estimated length of the queue may be 
25 necessarymordertodeterminemeappropriatewindowlengmtobeusedformequeu^ 

estimate, and to correctly estimate the length of queue. Such readjustment may be performed 
by changing the estimated length of the queue based on compensating for trends in the 
average arrival rate and for irregular departure rate. 

Another way to take into account this adjustment and readjustment, e.g., for a 
30 changing queue length, is to reduce (or extend) the length for each map by an amount equal 
to the length that the map is expected to be different from the average length or from the 
length of the last query. The longest of these adjusted lengths for the different maps is then 
preferably used as the estimate of the length. The estimated length may then be further 
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adjusted based on the trend so that it is referenced to the situation at the end of the window 
(last map) or to the average length during the window. 

Where me clearance is lmown o^ 
be further refined to remove the effects of variations in the clearance in particular cycles on 
5 the Poisson distribution. When using per cycle clearance information, the position of the 
farthest reporting probe is adjusted not only for changes in the expected length of the queue 
based on average differences as described in the previous paragraph but also for the known 
(and non-statistical) changes in queue length caused by variations in the number of vehicles 

exiting the queue for each cycle. 
10 In a preferred embodiment of the invention, the difference between the arrival and 

departure rates can be monitored. In one preferred embodiment, the rate is monitored by 
utilizing trend analysis of the lengths of the estimated queue length. Other methods. such as 
match filter analysis may also be possible. 

In the foregoing discussion, the map cycle time is preferably determined based on a 
15 predetermined timing, such as the timing of a traffic light, to give an optimum estimate. It 
should be understood however, that where the congestion is based on an accident or a merge 
or the like, the cycle time is not fixed andmay be chosen to achieve an acceptable system time 
resolution and/or length expected accuracy. The cycle time determines the motion per cycle 
(arrival and departure) and hence me resolution, wim higher rates reducing me a^^^ 
20 for congestion at a traffic light, the cycle time may be set at two changes of the light, with 
decreased performance, but utilizing reduced communication resources. 

The vehicle mapping and traffic reporting system of preferred embodiments of the 
present invention, provides wide latitude for distribution of computational resources within 
the system, between the individual vehicles, either when they are functioning as mapping 
25 vehicles or as receivers of information, and a central station. This distribution is possible, at 
least in part due to the simplicity of the computations required by the invention. 

In preferred embodiments of the invention, vehicles determine their positions either 
absolutely using two dimensional mapping or with respect to the focus of a congestion on a 
road using one dimensional mapping. As an end result, a map of areas of congestion, based 
30 on the responses to the queries is provided to the participating responding vehicles and/or to 
omervemclesandyorto Computerized 

in61ude an estimated time necessary to traverse the congestion. In addition to construction of 
the maps themselves this requires a number of relatively simple steps (as for example 
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described in 1 -5 above) and collation of the information. Depending on cost and bandwidth 
considerations, these computations may be distributed among (a) computers in the vehicles 
wMcharetrarisrmttmg^ 

a focus in the last few mapping cycles or a number of positions from the last few cycles), (b) 
5 in a central station which receives the positions and broadcasts them to all the vehicles in the 

pool and (c) in the receiving vehicles. 

Figs745illustrateasystemformedeteniunationandmappingofareasofcongestion 

in accordance with a preferred embodiment of the invention. 

Fig. 7 shows a response from vehicles in a first mapping stage, for identifying 
10 congested areas in which vehicles that are stopped or traveling below some speed have 
reported their positions preferably as described above, which are mapped as black rectangles 
in Fig. 7. A subsequent step can use one of two alternative methods, in accordance with a 
preferred embodiment of the invention. If a focus and branches can be identified according 
to the results of the first mapping step shown in Fig. 7, then a subsequent step is as shown in 
15 Fig. 9. If the focus cannot be defined the next step to be performed is described using Fig. 8. 

Fig. 8 shows a second, usually intermediate, mapping step in which vehicles (again 
mapped as black rectangles) in the vicinity of a focus of congestion (the circled area of Fig. 
7) have been asked to report their positions, preferably in slots, at a higher resolution (50 m) 
than in Fig. 7. Note that 50 m is occupied by about 10 cars. The gray slots do not refer to 
20 responses, but to areas of the jam. After determining that there appears to be a congestion in 
the area circled in Fig. 7, in a preferred embodiment of the invention, the system switches 
fromareabasedmappingtoafocusbasedm^^ 

described in Fig. 9. Note that up to this point there appear to be two vehicles reporting on 
branch 1, two reporting on branch 2 and one reporting on branch 4. However, Fig. 8 also 
25 shows the true length of the congestion which is under-reported in branches 1 and 2 and 

completely unreported in branch 3. 

It should be understood, in this and the following Figs., the vehicles are not identified, 
only their position, with respect to the focus, is mapped. Thus, when a vehicle position is 
identified in several successive maps, this information is presented for the information of the 
30 reader, but the ID's are not known to the mapping system. 

Fig. 9 shows a mapping response based on a congestion focus at the meeting of 
branches 1, 2, 3 and 4. Since the crossing of these branches is an intersection, it is assumed 
to be the focus of the congestion. In the response slots of Fig. 9, each of the branches (and not 
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the entire area) is mapped with a resolution of 5 meters (about one car length, except for 
branch 2 that is mapped, in this example, with lower resolution). In addition, where there are 
multiple lanes in each direction, the mapping (in an optional feature) also determines which 
lane(Ll, L2, L3, L4) the vehicleis in, with preferably, the same distance from the focus being 
5 mapped for each branch. For some lanes, for example turning lanes, the total length of lane 
that is mapped is shorter than other lanes. When positioning resolution of GPS or dead 
reckoning resolution is not sufficient to determine the lane of the vehicle, then lanes can not 
be distinguished. However, the position of the vehicle in right or left turn lanes may be 
inferred from the turn indicators of the vehicle activated by the driver. Moreover, lanes may 
10 be combined to present a single stream of congestion. This is sometimes preferred since 
merged responses of this type reduces the number of slots required. 

As can be seen, at the higher resolution, the single indication of a vehicle "c" in Fig. 
8 has resolved into three vehicles. Note further that the same number of slots were used in the 
mappings of Figs 8 and 9. When the lane information is not transmitted, more frequent 
15 updates, lower bandwidth use (freeing the remaining bandwidth for mapping other foci of 
congestion) or higher resolution may be achieved. 

It should be noted that if a focus (and its one or more related branches) can be 
identified from a low resolution query (such as that of Fig. 7), then the zooming of Fig. 8 may 
not be required. 

20 Figs. 10 and 11 show reporting by vehicles at some later time, using the 

reporting/mapping systems of Figs. 8 and 9. This later time is, however, close enough to the 
earlier report, that the congestion is assumed not to have changed significantly. Again, the 
reported vehicles are shown as black rectangles and the true extent ofthe congestion is shown 
in gray. Note that vehicle b has left the congested area (by passing the focus) and vehicles f 

25 and g have entered it. 

If only the reported positions ofthe vehicles as shown in Figs. 10 and 1 1 were used 
to determine the extent ofthe congestion, the length ofthe congestion in branch 2 would be 
under-reported and the congestion in branch 4 would be poorly reported. 

However, as indicated above, the structure ofthe congestion may be considered to be 
30 "quasi-stationary" at least over the time of several mappings cycles (e.g., one to a few 
minutes). With this in mind, for each ofthe branches, one may assume that the length ofthe 
congestion in each branch has not changed appreciably between the two (or more) reporting 
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periods and only the movement of the few reporting vehicles causes the maps to change. 
Altematively,thevehiclesthemselvesmay store their own position information for a number 
of cycles and broadcast either their largest distance from a focus (while in the slowdown) 
during the last few cycles, or may broadcast in a number of slots representing these past (and 

5 present) positions. 

Fig. 12 shows a "virtual" map of congestion at the intersection, produced in 

accord wimmeremapp^^ 

in particular with items 2) and 3) above. A virtual map includes not only present positions of 
the vehicles, but the entire estimated length of the congestion. The map of Fig. 12 shows the 
10 "known" congested area in black and an as yet unknown length of congestion in gray. The 
congestion lengths h a vebeendeducedbyutilizingthevehiclepositionsreportedinFigs.9and 
11. Note that the congestion length is known to a much higher accuracy than in any of the 
individual maps of Figs. 8-11. Further remapping cycles are preferably utilized to further 
improve the accuracy of system. 
15 It will be understood preferred embodiments of the invention, as described herein 

provides, with fewer reporting vehicles, improved resolution and accuracy of the length of a 
congestion, improved validity of congestion and frequent updates. By correctly choosing the 
timebetweenremappings and the number of remappings used in forming amap, the accuracy 
of the mapping may be varied at the expense of time resolution. 
20 The system described above may utilize a central decision maker which receives 

information from vehicles, plans the routing for each vehicle and then broadcasts a route or 
route changes to the individual vehicles. This type of system has the advantage that the 
routing for each vehicle takes into account the routing for the other vehicles and the control 
center, in computing the routings, can balance the routings to cause minimum delays or other 

25 optima*^ 

individual vehicles of their individual corrected routes. 

A second approach for routing systems which has been suggested is to have each of 
the vehicles compute its own route, based on some information about the present status of 
traffic which it receives from a central transmitter. While such systems require only a limited 
30 bandwidth, the routes computed by the individual vehicles cannot take into account the future 
effects of the routes of other vehicles. In further preferred embodiments of the invention, the 
actual congestion maps are also produced in the vehicles which receive raw, or somewhat 
collated information transmitted by the vehicles in congested areas. 
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An alternative system 190 of this type, in accordance with a preferred embodiment of 
the invention is shown in Fig. 1 3 . In Fig. 1 3, a plurality of local area transceivers 200 receive 
information from vehicles in regions surrounding transceivers 200. This information is 
preferably transferred to a concentrator 204 which receives information from a number of 

5 transceivers 200 and relays the information to a central station 206. Central station 206 then 
rebroadcasts the information (either as raw information or as maps or utilizing any other 
suitable format) to all the vehicles. Central station206 can also be used to generate the queries 
as well and then to broadcast multiplexed data containing queries and constructed traffic 
maps. When the queries and results are received by all receivers, they can utilized the latest 

1 0 broadcast results together with the stored (previous) queries to produce accurate maps, by the 

methods described above. 

In a preferred embodiment of the invention, vehicles are queried to transmit 
information as to which of a select number of troublesome intersections (including 
intersections already congested) that they expect to enter and when they expect to enter them. 
15 This information is preferably transmitted (in accordance with a query) in slots assigned to 
the intersections and estimated time intervals of arrival. Since more than one vehicle may 
expectto enter an intersection during a time interval represented by a given slot, in a preferred 
embodiment of the invention, a plurality of slots are assigned to each time interval and the 
number of vehicles is estimated from the number and percentage of slots in which a signal is 
20 received,usmgstatisticales^ 

criteria, responding randomly in one of the plurality of slots assigned to a time interval. 

Additionally or alternatively, the future development of existing slowdowns can be 
estimated from the prior development of the slowdowns, the rate of change of the length of 
the slowdown and the average speed of the vehicles that are within the slowdown. Such 
25 information can be made available to the vehicles based on comparison of the development 
of slowdowns which are detected by the methods described above. Such a method helps to 
construct a time developing map of intersections sensitive to traffic jams (trouble spots). 

Based on the estimates of the numbers of arrivals and times of arrival of the vehicles 
at the trouble spots, statistical information on future expected traffic jams is generated by the 
30 central station. However, in order to update the vehicles with real time expected traffic jams, 
the system has to perform periodic checks on trouble spots and to update vehicles with 
validated predictions of traffic jams, so they will recalculate their individual routings. This 
estimate of times of arrival may be based on a query based system, as described above, in 
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which the queries request information on the expected arrival time of the vehicles at various 
intersections, including at least those which are already congested. 

This recalculation of routes, broadcast of times of arrival at trouble spots and 
estimations of future traffic jams and slowdowns provides an adaptive refreshed process that 
5 uses current and predicted information, that gives each vehicle the information required to 
makeadistributedroutecalculationsystemeffectivein avoiding future problems, without the 
huge bandwidth requirements of central calculation of routes for the vehicles. 

Awntinuousprocessthatupdates^ 
predicted traffic jams, could be used With a distributed Dynamic Route Guidance (DRG) 
10 systemwhichdynamically selects the preferred routes. Such a DRG process performed, in the 
vehicle, should preferably help to alleviate congested roads in a synchronized way wherem 
different vehicles synchronize their processes of DRG. Synchronization is preferred to avoid 
too many vehicles taking the same route, causing congestion there. With such a method an a 
posteriori correction of try and fail processes can be used. This means that when vehicles 
15 receive predicted information about numbers of vehicles that are expected to pass a road or 
intersection in a given time, some of them should choose a "less preferred" alternative. This 
alternative may increase their travel time or distance by some factor; however, the overall 
result may be that traffic jams do not result or are less severe. If, on a second check of the 
predicted traffic, the jam still occurs, even less preferred alternative (i.e., ones that take even 
20 longer) are used by some of the vehicles, at least until travel time is equalized among the 
vehicles. 

Figs. 14 and 15 show two systems for connecting the mapping information, into car 
navigation systems (CNS) in accordance with preferred embodiments of the invention. 

Fig. 14 shows a simple system 220 in a vehicle in which an intelligent mapping 
25 transceiver 222 receives queries and slot allocations and sends position and/or other 
intaationinrespe^ 

224 via a standard interface, such as Japan's VICS interface, or the European RDS standard 
or other information interface formats. 

Fig. 15 shows a more sophisticated system in which the CNS may provide 

30 computationfaciliuesand^ 

positioning information (which the mapping transceiver uses to determine its position), GPS 
timing (which may be used as the master timing for the mapping system, with the slots being 
timed from a GPS timing signal which is common to all the vehicles and base stations), dead 
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reckoning positioning information (to improve the accuracy of the positions being reported) 
and/or map related information (so that the mapping system may provide traffic maps for the 
CNS). In addition, the mapping transceivers may receive and relay to the CNS traffic 
information from other sources which it may combine with its own traffic information before 
5 S endmgtomeCNS.Altematively,GPSinformationmaybedetermme^ 

receiver in the mapping transceiver (for example, 222 in Fig. 14), or from an external source 

(e.g. CNS provided data). 

Information may be sent by the control center to the vehicles to enable them to 

miiiimizeavemget^ 
10 consist of the above mentioned maps or maps of travel delay information at various 
mtersections-Thevehiclescanmenusethisiiifonnationtoop 

current and predicted maps may be used by the control center to send routing information to 
someofthevehiclesinorderto equalize traffic delays. In either event, the fest response of the 
mapping system allows for real time supervision, adjustment and continuous stabilization of 
15 traffic patterns with additional iterations. As described above, in a distributed system only 
prospective traffic patterns are broadcast by the control center and each vehicle calculates its 
own route. 

In many of the above embodiments of the invention, the system is triggered and/or 
synchronized according to a synch signal broadcast by the control station. Other sources of 
20 synchronization, which synchronize both the remote and control station, such as GPS 
received signals or other timing signals, can be used to trigger and/or synchronize the system. 

The invention has been described herein using examples in which the indication 
signals are transmitted in time, frequency or time and frequency slots. Other types of 
transmission slots are also useful in the invention such as frequency hopping and other spread- 
25 spectrum transmission slots. The term "transmission slots" or "slots" as used herein includes 
all these types of slots. In addition, while the invention has been described in a preferred 
embodiment thereof in which the positions of the probes are determined using the preferred 
query/ S lotresponsemethoddescribedabove,inotherpreferred embodiments of the invention, 
the actual reporting function may utilize other data transmission methods, such as Aloha, 
30 slotted Aloha or other methods known in the art. In such transmission methods, the distance 
from a focus, for example, is determined based on data specifying distance from the focus. 
Such memods may be useml when me percentage of 
more than about 5 percent. 
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The terms "comprise" "have" and "include" or their conjugates, where used herein, 
mean "including, but not necessarily limited to." 
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1 i. A method of estimating the position, in a ITS system, of the length of traffic 

2 congestion at a given focus, the method comprising: 
determining the positions of at least one vehicle farthest from the focus as a function 



3 



4 of time; and 



5 



estimating the length of the congestion based on the function. 



1 2. A method according to claim 1 wherein the position is estimated as the position of a 

2 vehicle farthest from the focus. 

1 3. A method according to claim 1 wherein the position is estimated as the position of a 

2 vehicle furthest from the focus during a given preceding time period. 

1 4. A method of improving the reliability of an ITS system comprising: 

2 determining the position of a plurality of vehicles; 

3 determining an indication of a traffic stoppage if more than one vehicle is stopped 

4 along a line of vehicles. 

1 5. A method of estimating the length of a road congestion, based on probe vehicles 

2 reporting characteristic values of their position to a receiver of a mapping system which 

3 processes the reports, the method comprising: 

4 (a) constructing a predetermined number of mapping samples, 

5 (b) determining in each mapping sample a position that relates to the position of 

6 a probe relatively far from a mapping focus, preferably a position close to the farthest probe 

7 position; and 

8 (c) selecting from the positions determined in step (b) the position which is the 

9 farthest from the mapping focus, thereby to determine an indication of the length of the road 

10 congestion. 

1 6. A method according to claim 5 wherein the position determined in step (b) is the 

2 position of the farthest probe from the mapping focus. 



SUBSTITUTE SHEET (RULE 26) 



WO 00/54143 



-44- 



PCT/mOO/00239 



1 7 . A method according to claim 5 or claim 6 wherein, after construction of a mapping 

2 sample, a response is transmitted to the reporters that disables transmitters that did not 

3 transmitareportwimmapreselectedrangemtheconstmctedmappmg 

4 disabled transmitters from continuing to report. 

1 8. A method according to claim 7 wherein the selected range comprises the position of 

2 the farthest probe. 

1 9 . A method of determining traffic motion and length of road congestion, in a system 

2 wherein probes, in response to a predetermined protocol, report characteristic values of 

3 their position to a receiver of the mapping system which processes the reports, the method 

4 comprising: 

5 (a) constructing a mapping sample that includes at least one of said reports; 

6 (b) selecting a range of said position characteristic values in which the farthest 

7 reporter from a mapping focus is identified in a mapping sample constructed in (a); 

8 (c) transmitting to reporters a response that according to a predetermined 

9 procedure disables transmitters that did not transmit a report within the selected range 
10 from continuing to report; 

! ! (d) receiving further reports and constructing a subsequent mapping sample; 

12 (e) repeating steps (a) to (d) according to a predetermined procedure; 

13 (f) selecting from the ranges selected in step (b) the farthest selected range to 

14 be indicative of the length of the road congestion; and 

15 (g) determining motion length toward a mapping focus by calculating a range 

16 characteristic value for a range in a mapping sample, subsequent to the first mapping 

17 sample, which includes the position characteristic value indicative of the closest position 

18 to the mapping focus and calculating the difference between the said range characteristic 

19 value and the range characteristic value of a corresponding selected range in an earlier 

20 mapping sample. 

1 10. A memodofcreatmgconditionswhichenableassessmentoftrafficmotionrate 

2 toward a mapping focus in a congested road and of the road congestion length at a certain 

3 time, wherein according to a predetermined protocol probes report characteristic values of 
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1 their position to a receiver of a mapping system which processes the reports, the method 

2 comprising: 

3 (a) constructing a first mapping sample that includes at least one of said 

4 reports; 

5 (b) determining a range of said position characteristic values in which at least 

6 one of said reports was identified in the first mapping sample; and 

7 ( C ) transmitting to reporters a response that according to a predetermined 

8 procedure disables transmitters that did not transmit a report within the selected range of 

9 the first mapping sample from continuing to report. 
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